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Abstract 

The Bethe-Salpeter equation for Wick-Cutkosky like models is solved in dressed ladder approxi- 
mation. The bare vertex truncation of the Dyson- Schwinger equations for propagators is combined 
with the dressed ladder Bethe-Salpeter equation for the scalar S-wave bound state amplitudes. 
With the help of spectral representation the results are obtained directly in Minkowski space. 
We give a new analytic formula for the resulting equation simplifying the numerical treatment. 
The bare ladder approximation of Bethe-Salpeter equation is compared with the one with dressed 
ladder. The elastic electromagnetic form factors is calculated within the relativistic impulse ap- 
proximation. 
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I. INTRODUCTION 



In the Quantum Field Theory the two body bound state is described by the three- 
point bound state vertex function or, equivalently, by Bethe-Salpeter (BS) amphtudes, both 
of them are solutions of the corresponding (see Fig. |l|) covariant four-dimensional Bethe- 
Salpeter equations (BSE)|1]]. Up to now most of the studies were restricted to the case when 
irreducible interaction kernel is approximated by (sum of) single particle exchanges. In this 
so called ladder approximation the scattering matrix is given by the sum of the generated 
ladders. It is known, that such approximation is not sufficient when more realistic models 
are considered 0, 0, To move beyond this approximation one is in practice confined to 
the use of some phenomenological Ansatze. (In hadronic physic, these Ansatze are very often 
made already at the level of two point correlators. For modeling of the gluon propagator in 
the context of BSE and Dyson- Schwinger equations (SDE), see for instance |^). 

Here we are considering some simple scalar models. An extensive review of BS studies in 
scalar theories (with at most cubic and non-derivative interaction) can be found in Ref . , 
0. The various improvements of simple ladder kernel have been considered, in particular, 
including the self-energy effects P] , or contributions from crossed box diagrams ITOf . The 
study of the influence on the bound state spectrum following from the infinite resummation 
of certain ladder and crossed-box diagrams can be found in the paper |Tl|. Furthermore, 



there is a number of interesting papers on the solution of Wick-Cutkosky models (with zero 
mass of the exchanged particle). These solutions employ various effective techniques like the 
point form of relativistic quantum mechanics [|12|, variational calculations [0 or the light 



front dynamics 

The standard approach to determine the spectrum and the BS vertex makes use of par- 
tial wave decomposition which reduces four-dimensional integral equation into the two- 
dimensional one. The alternative more recently exploited treatment is based on the 0(4) 
hyperspherical expansion . In this approach the BSE is transformed into an infinite set of 
one-dimensional integral equations. Notable advantage of this approach is a good numerical 
convergence and easy identification of excited bound states spectra. 

Very often, the ladder BSE's are solved with the help of the so-called Wick rotation 
|T^. However, the backward analytical continuation is quite difficult even for the ladder 



approximation, while for more complicated cases its proper implementation is unclear or at 
least highly non-trivial. 

In this work we follow the method of solving the BSE directly in Minkowski space H], |]T7 



in which the problems associated with Wick rotation do not arise. The method is based 
on utilization of generalized spectral representation for n-point Green functions in quantum 



field theory |jT9[. In this treatment the BSE written in momentum space is converted into 
a real integral equation for a real weight function with number of independent variables 
dependent on details of the model. We extend the earlier work |P first to the case with 
unequal masses of constituents. This then allows us to treat the ladder BSE in which 
all propagators (of constituents and of the exchanged particle) are fully dressed. This is 
achieved by the implementation of Lehmann representation of the propagator: 

G(p2) = J duja{uj)D{p;uj) , 

Dip;u;) = ^ ^ ; a{u) = R6{u - ml^iJ + aiu) . (1.1) 

p'^ — uj + te 

where (j{uj), which is a smooth function, nonzero above a threshold, is determined by the 
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Dyson-Schwinger equations (DSE). In the pole term we can take R = 1, a. choice corre- 
sponding to the conventional on-shell renormalization scheme {G has a unit residuum when 
momentum approaches a simple pole at physical mass). 

To account for the effect of self-energy we transform the momentum BSE to the form 
suitable for complementary solution together with the appropriate DSE for propagators. 
Note here, that the perturbative one loop contribution has been already considered in |^ 
and certain Euclidian version of this problem has also been investigated [^]. In qualitative 
agreement with |^ we have found that the critical value of coupling gives the domain of 
applicability of BSE (at least in its ladder approximation). The couplings below critical one 
allow only solutions for relatively weakly bound states. It is even more interesting that the 
effect of the propagator dressing on bound state spectra is rather small. In comparison with 
the bare ladder approximation the same binding energy is then achieved with the coupling 
smaller by about several per cents, even for values of the coupling close to the critical one. 

Clearly, when we take some or even all particle propagators dressed, the number of 
spectral integrations increases. Note here, that up to the rather exotic case of massless 
Wick-Cutkosky model the appropriate solution is not known analytically but must be found 
numerically. Mainly due to this reason we re-formulate the equation obtained by Kusaka et 
al P] and we offer the solution where the appropriate integral kernel is free of any additional 
numerical integration (see Appendix of Ref . ^ for the original solution) . The elimination of 
this numerical integration then not only improves numerical accuracy but also reasonably 
decreases the CPU time. 

To see explicitly the effect of radiative corrections we compare the derssed BSE results 
with its bare ladder approximation. We set the parameters of our model to that used in 
Refs. []T3| and [0 , to compare the bare ladder solutions to those obtained before in , |]T3 



Having solved equations for spectral functions, one can determine the BS amplitudes in 
an arbitrary reference frame. This makes this technique suitable for calculations of response 
to an external fields. In section IV. we briefly introduce the formulas defining the elastic 
charge form factor G{Q'^) in relativistic impulse approximation (RIA). Although the elastic 
form factor represents a simple dynamical observable, its Minkowskian calculation represents 
nontrivial task. For this purpose we consider the (massive) Wick-Cutkosky model given by 
Lagrangian gauged as follows: 




V{<P,) = [mi + ^703) 0^01 + + f <^3j 02 + ^^303 , (1-2) 

where covariant derivative is = d^ — ieA^. In our form factor calculations the effects of 
scalar dressing was not taken into account, since it would significantly increase computational 
complexity of the problem. Furthermore, it is assumed that e « g which implies that the 
interaction of the charged particle field 0i with the electromagnetic field can be treated 
perturbatively. As in Ref. 0, we have chosen the same coupling constant for interaction of 
the field 03 with the fields 0i and 02. The form factors were calculated for several sets of 
masses of constituent and exchanged particles. 
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II. DRESSED LADDER BETHE-SALPETER EQUATION 



The BS amplitude for bound state {4>i,4>2) in momentum space is defined through the 
Fourier transform of 

(O|T0i(xi)02(x2)|P) = e-^^-^(O|r0i(r72x)02(-r/ix)|P) =e-^-^ / -^^e~'P-^p, P) , 



(2.1) 

where X = rjiXi + rj2X2 and x = Xi — X2, so that xi = X + ri2X, X2 = X — rjix. Here 
Pi^2 are the four-momenta of particles corresponding to the fields 0i,2 that constitute the 
bound state (0i,02)- The total and relative momenta are then given as P = pi + p2 and 
p = {f]2Pi — V1P2), respectively, and P^ = M^, where M is mass of the bound state. Finally, 
P ■ X +p- X = pi-Xi +P2 ■X2- From now on we will put rji = ri2 = 1/2, which corresponds to 
the usual separation of center of mass motion for equal mass case, but can be also employed 
for unequal masses (although X is then not the coordinate of the center of mass). 

Introducing the BS vertex function F = iG{^G2^^, the homogeneous BSE for S-wave 
bound state reads 

r(p, P) = tJ ^ V{p, k; P)G,{k + P/2)G2{-k + P/2)V{k, P) . (2.2) 

The bound states appear as poles of the scattering matrix. The normalization condition for 
the BS vertex function follows from the requirement that the pole appropriate to a given 
bound state is a simple one: 



{2'nY5\p~k) (^GM)G2{P2) 



+iG,{p[)G2{p'2) (^^(P'^'^)) G,{p{)G2{p2) 



F(p,P), (2.3) 



where V{p,P) is the conjugate of V{p,P) and pi^2 = ±P + -P/2,p'i2 = ±^ + -P/2 (for the 
details see e.g. Ref. P). 

In this work we do not solve the BSE with most general irreducible scattering kernel 
V (the most general structure of the kernel V written in terms of its Perturbation Theory 



Integral Representation (PTIR) can be found in Ref. |§| or [19|). Here we restrict ourselves 
to the case of dressed ladder approximation with 03-exchange, in which V{p, k; P) = g^G^it), 
where t denotes the usual Mandelstam variable t = {p — k)"^. Note that for the bound state 
of particles (0i, ^2) only the t-channel interaction above is effective, whereas for the bound 
state of {(pi, 0j), i = 1, 2 one has to consider also possible u and s channel diagrams. In the 
present work we study the case (0i,02) due to its simplicity and leave the other cases for 
discussion elsewhere. Let us also recall that in the (dressed) ladder approximation, defined 
by the exchange of chargeless scalar 03, the photon coupling to 0i alone (one particle current, 
or in other words RIA) is by itself gauge invariant, when taken between the corresponding 
solutions of BSE. And the normalization condition (|2.3| ) in this approximation reduces to 
the condition G{Q'^) = 1. We will solve the BSE for massive constituents {mi 2 > 0) and 
ms > 0, the invariant mass of the bound state satisfies < P^ < (mi + m2)^. 

Taking the dressed kernel and the full propagators of constituent particles into account, 
the rhs of BSE (12.21) can be written as 



tg^ n / dai a(«i) / ■^D{k + P/2; ai)D{-k + P/2; a2)D{k - p; «3)F(A;, P) . (2.4) 
i=i'' \^^) 
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The interesting unequal-mass ladder case of Ref . pO[ is also described by Eq. ( p^) , although 
^1 = ''72 = 1/2 and p and k are not relative momenta. Since the dependence on momenta in 
( p.4| ) is explicit, one might always re-scale T{p, P) to proper relative momentum. 



The integral representation for the BS vertex function may be written as |T9| 

r(p,P)=fdz I da- — ^ . ^. (2.5) 

' J J [a- {p^ + zp- P + P^/A) -itY 

The positive integer n represents a free parameter without clear physical meaning. One 
can take advantage of this freedom of choice to pick up n so that the numerical solutions 
of integral equations for spectral functions are made more stable. The spectral functions 
p["'](a, z) for different n can be related by integration over a by parts. Kusaka et al p] choose 
n = 2 for their numerical solution of the BSE, we adopt the same value in this paper. 

The bare (symmetric) Wick-Cutkosky model corresponds to the choice: ai = a2 = rn^, 
the exchanged boson is massless (0:3 = 0) and no radiative corrections are considered. This 
model is particularly interesting because it is the only example of the nontrivial BSE which 
is solvable exactly [|l6l. For this model, there is no freedom in choice of n and (for the S-wave 



bound state) the expression (|2.5| ) reduces to the one dimensional PTIR: 

— [p^ + zp ■ P + _r^/4) — te 

Using a technique similar to the one used in ref. [§], the BSE can be converted to the 
following real integral equation for the real spectral function: 



""^a' , z') = \ dz / c/a V^["] («',/;«, z)p["] (a, z), (2.7) 



1 

P 

where we denoted A = g"^ / {A-k)"^ . The derivation is presented in the Appendix A, where the 
explicit expressions for particular choices n = 1, 2 are given. The central results of this work 
are expressions obtained for V^"'\ which are simpler that the ones presented in Ref. No 
additional integration is required which decreases the computer time necessary for numerical 
calculation. Besides, our formulas hold also for unequal masses of the constituents. The 
extension (together with above mentioned simplification) to the case of more complicated 
scattering kernel is not so straightforward, but we believe that it is possible. Note also, that 
due to the property of solid harmonic with respect to the integration over the momentum 
the presented procedure can easily be generalized for the bound state with non-zero spin 
(here, the total orbital momentum) pj. 

A bound state with equal composite masses is described by the vertex function F which 
is symmetric under the transformation P ■ p —P ■ p. In terms of the weight function this 
symmetry reads p{a,z) = p{a,—z). However, there are solutions that do not respect this 
symmetry even in the case of equal masses. These are usually called ghost solutions and 
the appropriate amplitudes have a negative norm. Such solutions are often considered to 
be nonphysical and it is supposed that they point at inner inconsistency in description of 
relativistic bound states within the BS formalism, at least in the ladder approximation. At 
this place, it is important to mention that the Lagrangian (|L^) describes the models that 
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are a subset of theories with potentials unbounded from below and in very strict sense they 
are discarded due to the vacuum instability. On the other hand one can assume, at least 
for sufficiently small couplings the existence of local minima of the potentials is sufficient to 
support of the existence of ground state of the theory. While in the large coupling regime, 
say for g/m » 1 no reasonable physics can be learned from the perturbation theory and/or 
from equation like ladder BS by itself. To conclude, we note that such ideas are supported by 
at least two facts. The ghost BS solutions do appear only for a large value of A. Furthermore, 
from the Dyson-Schwinger study we know that the scalar theory studied here makes sense 
only up to the certain critical coupling, see e.g. 0, [pifj. 

The important question arises, what is the validity of the full theory when renormalization 
is properly taken into account. Although the quantitative answer lies beyond the approxima- 
tions used in this article and requires more careful investigation, we make a simple attempt 
to find the domain in which self-consistent solutions of the BSE and Dyson-Schwinger equa- 
tions (within the framework of reasonable approximations) exist. Furthermore, we inquire 
an influence of scalar propagator dressing on the solution of the BSE for the bound states. 



III. DRESSING PROPAGATORS BY THE DYSON-SCHWINGER EQUATIONS 

The solution of the DSE for the scalar models with the help of spectral decomposition will 



be discussed in detail in our forthcoming paper [^]. Here we give only brief presentation of 
the DSE in bare vertex approximation, their renormalization, re-arrangement in terms of the 
spectral function and some properties of solutions, important for our further discussion of 
the BSE. In this section we assume < mi + m2, so that the propagator of the exchanged 
particle $3 has an isolated physical pole. 

By dressing of the scalar propagators in our study of BSE we mean only the dressing 
due to the "strong" interaction between scalars, the coupling to the electromagnetic fleld is 
neglected. Let us now write the strong interaction part of our Lagrangian (|1.2|) in terms of 
bare, unrenormalized quantities (fields and coupling constants), labeled by subscript "0": 

^strong = —901 </'oi</'O10O3 — ^ 002003 • (3.1) 

In the previous section we have chosen the strength of of both couplings to be the same. 
Here, we distinguish the bare couplings, anticipating that they are renormalized by different 
amounts (see below). 

The kinetic terms are parametrized by the unrenormalized masses moj. These masses 
undergo the infinite mass renormalization 

fT^li = "^i - , i = 1,2,3. (3.2) 

To re-scale the residuum of the full propagators to unity, we will complement the infinite 
mass renormalization by the finite (since the model is superrenormalizable) renormalization 
of the fields and coupling constants 

4>oi = \/Zi4>i, ^ = 1,2,3, goi = —^j=gi, i = 1,2 . (3.3) 



That is, we will employ below the on shell renormalization scheme in which the propagators 
have unit residua when momentum approaches its mass shell value — > m^. 
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In this paper we consider the Dyson-Schwinger equations in the simplest approximation 
in which the proper vertices are replaced by the bare ones Toi = Qoi- Then, the DSE in their 
unrenormalized form read 

Go^' (P) = P'- ml - no.(p2) , ^ = 1, 2, 3 

noi(p^) = igli 7—4 Go3 (p - q) Goi (q) , i = l,2 

J (27r) 

no3(p^) ='J-0f {P - q) Go. (q) , (3.4) 



1=1,2 



where Go{p) is the Fourier transform of the full unrenormalized propagator Got (2; ~ u) = 
< 0\T(f)oi{x)(j)Qi{y)\ > and Iloi is the corresponding self-energy. 

Under the field strength renormalization the propagators scale like Goi = ZiGi. Multi- 
plying the equations for G^l in (|3^ ), defining Ilj = ZJloi and making use of (|3.3| ), one gets 
the re-scaled DSE 

Gr\p') = Z,(/-m^J-n,(/), 
n.(/) =^9^|^,Gsip-q)Gdq), ^=1,2 

Mp') = EG.ip- q) G, (q) . (3.5) 

(2vr) i=i,2 

The renormalization of proper self-energies proceeds by double subtraction: 



U,^{p') = n,(/) - n,(m^) - ip' - mf) . (3.6) 



Sm^^ = Ui{m^i)/Zi , Zi = l+ j ^ |p2^^2 , (3.7) 



Identifying the appropriate renormalization constants ( |3.2| , |3.3|) 

dp^ 

we can immediately write the full propagator in terms of finite physical quantities 

GT' (P) = / - ml - n,^(p2) , ^ = 1, 2, 3 . (3.8) 



The DSE for the renormalized self-energies are given by ( p.5|) and subtraction (|3.6| ). 

For the purpose of our BS calculation, we now fix the renormalized couplings and masses 
as follows 

m 

91 = 92 = 9', mi=m2=m; 1713 = — (3.9) 



where g is coupling constant from (|1.2|) . That is, we will compare the solutions of the BSE 
for the bare and dressed ladder kernel taken for the same numerical value of unrenormalized 
and renormalized coupling constant, respectively. The masses are fixed to allow comparison 
with some of the results of Refs. H, |T^. 

Now, it is straightforward task to evaluate the spectral representation of the renormalized 
self-energy. Lehmann representation (with unit residuum) for Gi reads: 

G.{p')= Hds ^ ""^'^ . , ais) = S{m^ - s) + a{s) , (3.10) 
JO p'^ — s + te 



Notice that functions ai,i = 1,2, have the thresholds at rrii^th = {f^ii + '^s)^ = 2.25m^, 
whereas the function cxs has the threshold at m^^th = (^i + ^2)^ = 4m^. Analogously, for 
the for self-energies: 

n..V)^r /"'°>. f -'";>: . (3.11) 

The spectral representation for explicitly satisfies nij(m^) = n'^(m^) = following from 
( p.6|) . Rewriting now the relation between G and 11 in the form G = D + DUG {D being 
the free propagator with the physical mass) and taking its imaginary part, we arrive to the 
first relation between the spectral functions a and p: 

(Ti{uj)p^X(^) , ai{a)p^X^) 



(ul) — r J uo — a 



1,2,3 



(3.12) 

where P J stands for principal value integration. All the functions in ( p. 121 ) are positive and 



regular above the perturbative thresholds and identically equal to zero elsewhere. 

Substituting the spectral representations (|3.10|) into the DSEs (|0| ), making the sub- 
traction as in ( p.6|) and comparing to the Ihs in the form of ( 3.11|) , one gets after lengthy 
algebra: 



p^X^) = ^ J dadl3B{a,P;uj)as{a)cri{/3) , z = l,2, 

P-kX^) = ^ ^2 dad(3 B{a, j3]Lo)ai{a)ai{l3) , (3.13) 

1=1,2 

where A = g'^ / (47r)^ and the function B{a, (3; u) is related to the Kallen function A as follows: 

X{a,p,u) = + + LU^ - 2ap - 2auj - 2puj, . (3.14) 

Before the numerical treatment the explicit integration- separating the 5-function parts of 
Lehmann weights a- has to be performed. 

Equations ( |3.12y3.13| ) constitute the closed system of integral equations for spectral func- 
tions which can be solved numerically by iterations without any additional approximation. 
So obtained dressed propagators have been used when solving the Bethe-Salpeter equa- 
tion. The results are discussed in the section V. Before leaving this section we review some 
important features of our solutions of DSE. 

The behavior of the imaginary parts of propagators- the Lehmann functions ai{a)- for 
fields $1,2,3 is shown at Fig. ^ 

The renormalization constant Zi are calculated from the relation 

Z, = l- fda . . (3.15) 

J [a — mfY 

From the Fig.^ we can see that the field renormalization constant Zi,2 changes sign from 
positive to negative at some critical point \crit = QlritI {^'^'''^Y — 1-5 ± 0.1, where the error 
reflects the difficulty of making the numerical estimate of the value for which the solution 
cannot be found and the dimensionless coupling is defined as A = \/m? . We did not find any 
numerical solutions of DSE for couplings larger than Xcrit- It is reasonable to suppose that 
the quanta associated with the fields 0i,2 do not describe physical particles when A > Xcrit- 
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IV. ELASTIC ELECTROMAGNETIC FORM FACTOR 



The electromagnetic form factors parameterize the response of bound systems to external 
electromagnetic field. The calculation of these observables within the BS framework proceeds 
along the Mandelstam's formalism |p2|. For the elastic scattering on the S-wave bound state, 

f 

element in terms of the single real form factor G{Q'^) 



[Pf = Pf = M^) the current conservation implies the parameterization of the current matrix 



G^^{Pf,P,) = G{Q'){P, + Pfr- (4.1) 

The elastic electromagnetic form factor G{Q'^) depends only on the square of photon incom- 
ing momentum q and we use the usual SLAC convention Q"^ = — g^, so that is positive 
for elastic kinematics. 

The matrix element of the current in relativistic impulse approximation (RIA) is dia- 
grammatically depicted in Fig. ^. In this paper we are not taking into account the dressing 
of the scalar propagators when calculating the charge form factor. Then, the matrix element 
is given in terms of the BS vertex functions as 



G^{P + q,P) = I I ^f(fc+|,P + g) 



DiPf, ml)j>^ipf, p,)Dip,; ml) D{-k + P/2; ml)] T{k, P) , (4.2) 

where we denote P = Pi and jf represents one-body current for particle 0i, which for the 
bare particle reads ji{pf, Pi) = +pf , where PiiPj is initial and final momentum of charged 
particle inside the loop in Fig. Q i.e., Pi = k + P/2,pf = q + k + P/2. 

We have already mentioned in the previous section, that if the vertex functions F are 
solution of the BSE with a kernel corresponding to exchange of single chargeless particle, 
the RIA defined above is by itself gauge invariant and the normalization condition for the 
BS amplitudes is equivalent to the normalization G(0) = 1. 

The main result of this paper, as far as charge form factor is concerned, is re-writing 
of the rhs of Eq. (^4.2|) directly in terms of the spectral weights of the bound state vertex 
function. It allows the evaluation of the form factor by calculating the integral of nonsingular 
expression, without having to reconstruct the vertex functions F(p, P) from their spectral 
representation. The derivation of this integral involves some lengthy algebra and is relegated 
to Appendix B. 



V. NUMERICAL RESULTS 
A. Bare Ladder BSE 

We have solved the bare ladder BSE for symmetric (mi = m2 = m) scalar theory with 
bare ladder kernel 

V{p, k, P) = Vip -k) = 3 , (5.1) 

[p — /e)^ — m2 

and bare constituent propagators Gi{pi) = Di{pi,mi) by iterations of the integral equation 
for spectral functions. The standard procedure was followed: after fixing the bound state 
mass (P^) we looked for the solution by iterating spectral function for fixed dimensionless 
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"coupling strength" A = (7^/((47r )^m^). If the iterations failed- measure being both the 
difference of the rhs and Ihs of the integral equation and deviation of the auxiliary normal- 
ization integral from pre-defined value- we were changing A (halving intervals of successive 
guesses) till the solution was found. 

In case of Wick-Cutkosky model the one dimensional integral equation ( |A24| ) was solved. 



Although the solution of this one- dimensional integral equation could be found by the in- 
version of its discretized form, we have tested the iteration procedure (used later also for 
massive exchange). The equation was discretized by numerical Gauss integration, it appear 
that it is sufficient to take 40 Gauss point (though the number cited in Table 1 are obtained 
with 98 points). It is known [^] that = corresponds to A = 2 from which our result 
slightly deviates in the fifth digit. We have also reproduced (up to four published digits) all 
results for Wick-Cutkosky model from |]TB[. In the Fig.^ the weight functions are plotted 



against spectral variable z for several fractions of binding t] = \/P'^/2m . 

For the massive scalar exchange the two-dimensional integral equations ( |A31| , [A35D were 



solved. We have found (in agreement with P) that numerical errors are about one order of 
magnitude bigger for n = 1, hence n = 2 is preferable and only the results with this choice 
are discussed below. 

For numerical solution we discretize integration variables a and z using Gauss-Legendre 
quadratures (with tangent mapping from < — 1,+1 >^< a^m, oo > for a ). The equation 
is solved on the net of = Nz*Na points which are spread on the rectangle (— 1, +1)* 



{amin, oo). The value amin is given by the support of the spectral function (see the appendix 
A). We have not optimized the grid during the iteration procedure as it was done in the 
study p. Instead, we have solved the equation for several different numbers of grid points 
while keeping fixed the ratio of Na/N^ and then extrapolated the results to the "ideal" case 
with Na = Nz = oo. Examples of numerical convergence for some cases of bound states are 
presented in the Table 2. In the Table 1 we compare our results for = m/2 with those 
ofref. §. 

Below we show the dependence of charge form factor on the parameters of the model: 
on the range of interaction characterized by the inverse mass of exchanged meson rris and 
on the strength of forces which bind the particles together. To calculate the form factors 
we first have to solve the BSE for chosen sets of parameters. We vary the parameters as 
follows: 

1. First we solve the BSE for several bound state masses keeping the ratio m^/m 
fixed 

2. Then we vary the mass of the exchanged meson, keeping the masses of all studied 
bound states fixed (our choice is ?7 = \/P^/{2m) = 0.5) and determining the corre- 
sponding coupling strengths A. 



Where the independent numbers were available |]T3|, they agree with our results (see Table 
3). If the mass of exchanged particle becomes small (but nonzero) the convergence of 
our numerical procedure becomes somewhat poorer and more sensitive to the initial guess. 
For illustration the weight functions p^^^ is plotted in figs.(PJ^) for two different values of 
exchanged mass. 
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B. Dressed Ladder BSE 



In this subsection we finally discuss numerical solutions of the BSE including the dressing 
of the propagators. We introduce the dressing by two steps, switching it on first only 
for the exchanged particle and in the second step also for constituents. We should point 
out that since the solution of the DSEs in the bare vertex approximation (by which we 
dress the propagators) breaks down for coupling constants larger that Xcrit = 1-5 we can 
consider to only rather weakly bound states: for the bare BSE A = 1.5 corresponds to 
= /p2/(2m) ~ 0.78. 

The propagator dressing of the exchanged particle in the one loop approximation was 
already considered in the Ref. IQ. We go beyond the one loop approximation and determine 
the continuum part of Lehmann weight a-^ from the DSE ( |3.12| ) with the same value of the 
coupling constant. That is, we insert into Eq. (^]^ ) the dressed kernel 

Gsip -q)= r du- (5.2) 

with the pole situated at m^, = m/2. As noted above, the constituent propagators are at 
this stage left undressed in the BSE, although they all self-energies have been taken fully 
into account in DSEs. The integration over uj ( |5.2| ) in the BSE kernel was performed using 
Gaussian quadrature with 16 points. Including the kernel self energy slightly decreases (by 
at most few per cent) the mass of the bound state, even for A ~ Xcrit.- 

Let us point out, that the kernel of Eq. ( p.7|) in the dressed ladder kernel approximation 
is free of any singularities. The accuracy of the numerical solution is comparable to the the 
bare ladder case. For example, we A = 0.734 for t] = 0.95 for the grid of 32 * 32 points and 
A = 0.749(0.752) for the grid of 64 * 64(96 * 96), the convergence is similar to to the case of 
bare ladder (see the Tab.l). The extrapolated (to very large grid) values of X's for fractional 
binding rj = 0.999, 0.99, 0.97, 0.95 are showed in Fig.|. 

In the next step we have included the self energies of the constituents. As we shall see the 
effect is relatively small for A << Xcrit, but increases rapidly as A — > 1. As in the previous 
case the Lehmann weights have been calculated from the DSEs solved for the same value 
of the coupling. As the first guess we have used the solution of the BSE linearized in cr{a), 
i.e. with only one propagator dressed. This guess is rather close to the exact solution for 
T] < 0.9. 

The constituent particle in a weakly bound system (r^ ~ 1) live near their mass shell. 
Therefore, one can naively assume that the values of coupling for such weakly bound state 
should not be strongly affected by dressing of constituent propagators. For deeper bound 
states we have found that the effect of the dressing of constituents is much larger than that 
due to the dressing of the kernel (exchanged particle) Fig.|[ The couplings for fully dressed 
BSE are not determined with the same high accuracy as those for the bare ladder BSE, 
since the grids are not optimalized for very different ratios of a's which appear in the kernel 
of the BSE. 

C. Charge form factor 

Various form factors are extensively studied in scalar theories like Wick-Cutkosky model 
(see for example [Q). In these studies the dependence of the form factor on the binding 
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and on the range of the "strong" interaction has been considered, therefore we perform a 
similar calculations in our formalism. 

In the approach adopted in this paper (employing the spectral representations in the 
Minkowski space), the bound states masses and corresponding vertex functions can be ob- 
tained with good accuracy and in reasonable CPU time. Unfortunately, the calculation of 
the scalar form factor as outlined in appendix B leads to more complicated results ( |BliyBl2D . 
Even if one would be able to perform analytically all additional Feynman integrations, the 
formula ( |B12| ) still involves the four- dimensional integration over the spectral variables. 
We are taking also the integrals over four Feynman variables numerically with the help of 
Gauss-Legendre quadrature, taking the number of points for each of them equal to the one 
for spectral variable z. Since relatively small number of integration points (from 16 to 40) 
was taken for each integration, the presented results have to be viewed rather as an estimate 
of form factor behavior. One can always refine the grids at the expense of longer CPU time. 

We have also compared our results to choose obtained in the Gross (spectator) formalism, 
choosing the "scalar deuteron" parameters (see @). In analogy with the real deuteron, the 
parameters are chosen as: m^ jm = 138/938.9,7/ = (2 * 938.9 - 2.3)/2 * 938.9 ~ 0.9988. 
The bound state vertex functions were found by solution of the Gross and BS equations. 
All phenomenological form factors introduced in have been "switched off' (the limit 
A's — » oo are taken in the "strong" form factors) when calculating the Gross wave function 
and the bound state current. The e.m form factors were calculated in the spectator RIA 
and as described in appendix B (using the grid 32^), respectively. 

The form factors for several bound states listed in the Tables 1 and 3 are presented in 
Figs.| and |T^, respectively. In Fig.|^ the ratio of the exchanged and constituent mass is kept 
constant rris/m = 0.5 and the mass of the bound state M = y/P^ is varied. In agreement 
with physical expectations one sees that as the bound state become more tight the elastic 
form factor increases. (For the infinitely bound point system is should be equal to unity, 
even our deepest bound state are of course still only approaching this value.) We put into 
the same plot for comparison also the scalar deuteron result. How the form factror changes 
with the mass of the exchanged particle is shown in Fig.^. We included several states for 
which M = m, bound by the one-boson-exchange potential of the range r = which 
is varied. Two other systems, first with rj = 0.4; ma = 0.25 and the scalar deuteron, are 
added for comparison. From both figure we can conclude that the behavior of form factor 
is determined by the strength of the interaction and its range rather differently for various 
M. The range of the interaction is more significant for a weaker ones. This agrees with 



conclusions of ref. [|10|, which have also compared our results for small exchanged mass 
7] = 0.8, ma = 0.15m with Wick-Cutkosky model prediction for 77 = 0.784, ma = (Tables 2 
and 1 of [|10[, respectively), and found only slight difference in the range = (0, lOOm^). 



VI. CONCLUDING REMARKS 

The spectral representation was employed for solving the Bethe-Salpeter equation in 
(3+1) Minkowski space. The new analytical formula for the integral equation kernel has 
been derived. 

The method is efficient solving both the bare and dressed ladder BSE. Solution of Dyson- 
Schwinger equations for propagators leads to appearance of critical value of the coupling 
constant, beyond which the solution collapses. This restricts substantially the region in 
which the effects of dressing can be studied. Since the coupling is rather weak, the dressing 
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leads only to moderate decrease of the bound state masses: even close to the critical value of 
the coupling the fractional binding of the bound state of the dressed BSE is smaller than the 
corresponding one for the bare BSE by at most 15 per cents. As an example of application 
of the obtained vertex functions, we calculated the elastic electromagnetic form factor. 

To further develop the method, it would be interesting to extend it to more complicated 
BS kernel: trying to include the cross boxed contributions, s and u channel interactions 
etc. It is already known that the "spectral" approach used here is suitable even for more 
complicated systems, for scalar QED see Ref. |TB[. One of our future goals is to manage the 
complication due to fermionic degrees of freedom. 
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APPENDIX A: KERNEL FUNCTIONS 



In this appendix the real integral equation for the BS vertex weight is derived in detail. 
The PTIR form for scalar bound-state vertex reads 

r(p.P)^/d. / <i"^0^ . (Ai) 

where p^^^{a, z) is the real PTIR weight function for the bound state vertex function, and n 
is a dummy parameter. The function F is given by |jl9 

z;p,P) = a - (p"^ + zp- P + P^/A) - ie = a - f{p, P, z) - ie . (A2) 

The support of p["](a, z) can be determined in general (see |]l9l) for arbitrary interaction. In 
our case of one-boson-exchange interaction kernel one gets a bit higher amin • We discuss 
our treatment of the lower bounds below. 

The following procedure is straightforward but a bit exhaustive. The rhs of the BSE ( p.2|) 
with the kernel given by the exchange of the single (dressed) particle $3 has to be re-written 
in the form allowing to extract the integral equation for the spectral function p["](a, z). The 
integrand contains the two "constituent" propagators, the denominator F{a, z, k, P) from 
the spectral representation of r(A;, P) and the propagator of the exchanged particle (all other 
factors will be skipped for a while for the sake of briefness). 

Using the Feynman parameterization technique we first write 

Dik + P/2; a,)Di-k^ P/2; a,) = 1 / |,,. _ ^) _ ^^j. . 

M^(,) = ^ + ^„. (A3) 
Next the denominator of (lAll) is added: 



1 1 



D{k + P/2;ai)D{-k + P/2;a2) T{n + 2) r, r ,^ (1 - i)^ 



dr] I dt- 



n-l 



[F{a,z;k,P)]- 2Tin) 7 '7 [R - f{k, P, z') - ie]-+^ ' 

R = at+{1- t)M^ , (A4) 

where z' = tz + {1 — 1)7]. Now, we include the propagator of the exchanged particle, the 
integral over d'^k and factors ig'^, defining: 

tonn. = ^9 j [F(a,z,p,P)r <^^' 



. 2 r(n + 3) 
21 [n] 

d^k 



11 1 

J dr]Jdt{l- t)e-^ j dxx''+^ Ik 



{27ry 



X 



k'^ + 2k-Q-{l- x)p' - ^P^ + (1 - x)a3 + xR-ie 
i T{n + l) 1 1 



-(n+3) 



(47r)2 r(n + 3) x"+i(l - a;)"+i [A - f{k, P, z') - ie]"+i 

R ct'^ X ^ , . , 

A=- + --7^ A6 

1 — X X [1 — X) 
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where Q = {1 — x)p — xz'P/2 and S = — f— -P^- Since z' lies in the interval < — 1,+1 >, 
< 5" < (mi + m2)^/4. Interchanging the integrals over rj and t with the help of: 



-1 



1± z' 

r± = and = + (1 - i)77 , 



l±z 

and introducing A = g'^/^in)^ we get 

1 1 



(A7) 



W^Af/d.'/^|^i:ew.-/))/ 

_1 5— It n 



n—l 



[F{A,z',p,P)]-+^' 



(A8) 



Let us separate the t dependence of F{A, z';p, P). First we substitute for rj into the definition 
of M^, which yields 



ai + a2 ai — a2 \ ai + 0:2 ai — 0:2 / 



+ 



z + 



+ 



z . 



Next, we introduce the notation (indicating explicitly the t-dependence of A and R): 



1 — X X 1 — X 



1 — X X 



R^t) = at+{l- t)M' = J{a, z)t + ^^4^ + ^^—^ z' 

^ Oil + 0:2 Oil — OL2 

J[a,z) = a z, 

in which the dependence of F{A{t), z',p, P) reads 

F{A{t), z';p, P) = ^-^^ t + F(A(0), z';p, P) . 
Since t-dependence of F{A{t), z';p, P) is hnear, the integral over t can be taken: 



/ 



and hence 



F{A{t), z';p, P)"+i nF(A(0), z';p, P) [F{A{t), z';p, P)]" 

Q{s{z-z'))T: 



Idddf - 2 7 i n - xY+^ h 



-1 

Using this result, the BSE can be written as follows: 



F{A{0),z';p,P) [F{A{T,), z';p, P)]- ' 



(A9) 

(AlO) 
(All) 
(A12) 

(A13) 

(A14) 

(A15) 



/ dz' / da' , , 

I, J, [F{a',z',p,P)] 



J a J da J dzp^''\a, z) Idddf ■ (A16) 
j (T = j daiai{ai) j da2a2{a2) j da^a^{a^) , (A17) 



0:3,1 
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where a{ai) are the Lehmann functions of the dressed propagators, reducing to the (5-function 
if the dressing is neglected. Obviously, the rhs is still not quite in the desired form, both 
F{A{0), z',p, P) and F{A{Ts), z',p, P)" are functions of momenta p and P. 

Below we re-write the kernel as a sum of several fractions, which after substitution a' = 
A{T),T = 0, T± would allow to use the uniqueness theorem and extract the BSE in the 



spectral form. This is possible only if the integrals over a' on both right and left hand sides 
of |A16| are taken over the same intervals. To show this it is first necessary to prove that 
i?(T) — 5* > 0, for T = 0,T±, since the functions A{T) —>■ +oo for a; — *• and x —>■ 1 they 
have on for < x < 1 the minimum equal to: 

ArmniT) = [^sj R{T) -3 + 771^^ + 3 = R{T) +ml + 2m3 {R{T) - S) . (A18) 

next, one has to show that these lower bounds are not in conflict with the lower bound for 
a' (and the same lower bound for a). It is simple for the undressed equal mass case, when 
the condition above taken for T = actually defines the lower bound for a in the following 
form ^ 

a > (^^Jm? - 3{z') + mg) + S{z') , (A19) 

and the same for a,z ^ a', z' . Since R{T±) depend also on a and z and since for the equal 
mass case R{T±) > R{0), the next two constrains clearly conform with the lower bound for 
a' and one can extract from them the upper bound for the integration over a: 



2 1 
a < m + — 



[^Ja' - 3{z') - mg) +3{z') 



2 



(A20) 



(compare to eq.(A6) of P|, where one bracket seems to be misplaced). For the unequal 
mass case or for the dressed propagators this analysis is much more complicated, mostly 
due to the fact that now for some combination of parameters it can occur R{T±) < R{0). 
The necessary condition R{T) — 3 > can again be proven in this case (though after much 
longer algebra), ensuring that the common lower bound for a's exists. But we could not 
resolve the conditions (|A18 ) analytically, they are treated numerically. 



Now, let so go back to eq. ( |A16|) and proceed by considering first the simple case of the 
symmetric Wick-Cutkosky model. 



1. The Wick-Cutkosky Model 

In this model the constituents have the same masses mi = m2 = m and the mass of 
the exchanged particle is zero (all propagator dressings are neglected). Then for the S- 
wave ground state vertex function, the spectral function depends only on variable z and the 
denominator enters with power ?7, = 1 (|27 



p^''\a,z) ~^5{a-m^)p{z), (A21) 

and in Idddf we replace 

n ^ 1 , ^ m^ R{t) rr? , J(a, z) ^ , 



VP? X 



A{t)^A=- 3., 

\ — X \ — X 
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Taking the integral over x we find (recall that for equal masses < S < m^): 

^^^^^ ^ 2 I m^ + f{p,P,z')-te ' ^^^^^ 

Comparing both sides of the BSE and using the uniqueness theorem [jl9|, we get the well- 
known (see e.g. 0) integral equation for pW(2): 

1 

p{z') = X [ dzV^^\z' , z) p{z) , (A23) 



VKf,' z) - T Q{s{z-z'))T, 



s=± 

Although its solution is known analytically even for exited states, the energy spectrum still 
has to be found numerically (up to the = corresponding to A = 2). For the purpose of 
numerical treatment, this equation is usually rewritten in the form: 

^ = Vo{z')- ldzp{z)V^iz',z) 
-1 

1 , ^e(z' - z) + ^e(z - z') 

2(m^ — S) 2{m'^ — S) 

where the temporary auxiliary normalization condition / dzp{z) = 1 was imposed (i.e., if 
used in further application, the vertex function would have to be re-normalized in accordance 
to (pp). 



2. The BSE forn = 1 

We will now bring eqs. ( |A15|JA16|) to the desired form for the particular choice n = 1. 
For the numerical solution this value is not the most suitable one. But since the formal 
manipulations are in this case the simplest, we treat it first for methodical reasons. 

For n = 1 the integrand of ( [A15| ) can be decomposed into the sum of simpler fractions 
with the help of 

Q{s{z-z'))Ts ^ 1-x ^ 1 /^24) 

^^F{A{0),z';p,P)F{A{T,),z';p,P) J{a,z) ^ F{A{T), z'; p, P) ' ^ ' 

where we have introduced a shorthand notation 

E f{T) = /(O) - 9{z - z')f{T^) - e{z' - z)f{T.) . (A25) 
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Notice that the Ihs of ( [A24| ) is nonsingular for J(a, z) = 0, so when this happens the rhs 
behaves hke 0/0, which calls for some caution in the numerics. So, the integral Idddf now 
reads 

A f dx 

^"""^t" = 1) = 2JRI) (l-.)FiAinz';p,P) ■ <^2«' 

In the last step we introduce the spectral variable a' = A(T) and use the dependence of 
A{T) on X to convert the integration over x into integral over a'. Picking up explicitly the 
x-dependence of A{T) we can write: 

1 — X X 



R{T) - S as 1 f . as 
1 — X r x"^ 1 — X V x^ 



5{a'-g{x)) = J2 \ ,} .. S{x-Xi), g{x±) = a' , 

i=± \9 [Xi)\ 

x±{T) = ~ ~ P ^ ^ , D = {R{T) -a' + a,f - Aa,{R{T) - S) . 

E(x±, S,a') = a' -^-S. (A27) 

x± 

With the help of these relations we get 

^ /■ dx r da' ^ ^ 9{xi{T))9{l - Xi{T))9{D) 

{l-x)F{A{T),z'-p,P)~ J Fia',z';p,P)^ \E{x,{T), S, a')\ 



Using this result in ([A16D , one gets from the uniqueness theorem the integral equation for 



the BSE structure function is obtained: 

= A [a [ dz [ daV^^\a,z,a',z')p^^\a,z), (A28) 



Notice, that x±(T) depend on a',z' (which are fixed from the Ihs), ai,a2,a3 (which are 
fixed when the dressing is neglected), P^ in S (which is given by the binding energy of the 
system) and for T 7^ also on a,z (through J{a,z) in R(T±)). This allows to recast the 
integral equation into the form more convenient for numerical treatment. 

a. BSE without propagator dressing (for n = 1) 

If the propagator dressing is omitted, mf and the corresponding BSE for the spectral 
function reads 

pW(a',2') = X JdzJ daV^^\a',z',a,z)p^^\a,z) , 
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z\ a, z) 



1 d{x,{T))B{\-Xi{T))e{D) 
2J{a,z) V.ti \E{x,{T),S,a')\ 



As mentioned above, for T = the kernel depends on a and z only through J(a, z), hence 
it is convenient to pick up this case from the sum over T and re-scaling 

p^^\a,z) = J{a, z) p^^^a, z) , 

V^'\a',z',a,z) = y V^'\a',z',a,z), 
J [a, z) 

and imposing the auxiliary normalization 

dz J apW(a,z) = 1, (A30) 
rewrite the BSE in the non-homogenous form: 

p[^l (a', z') = X Vi^^ («', z')-X JdzJaY. V}^^ (a', z' , a, z) (a, z) , (A31) 



2 J(a', z 



^^(x,(T))e(l-x,(T))e(D) 



|E(x±(T),5,a')l 



T=0 



2J(a',;z' 



^(x,(T))g(l-x,(T))g(D) 
|E(x±(T),^,a')l 



The ^ functions in the kernel impose the proper bounds on a and ensure that the rhs 
contributes only for a' from the support of p[^](Q;', z'). 



3. BSE forn = 2 



Now, we would first describe the necessary modifications for the choice n = 2 which was 
used in actual numerical calculations in this paper. Going back to ( |A15|) we can for this 
case write: 

e{s{z-z'))T^ 



E 



iF(A(0),^',p,P)F(A(T,),z',p,P)2 

1 — X ^ r T 

J {a, z) ^ ^ 



+ 



1 — X 



\F[A{T),z'-p,P)]^ J{a,z) F{A{T),z';p,P) 
Then, the integral Idddf can be cast into form: 



A 



IDFFF 



' dz' C E 



T 



+ 



1 — X 



1 



2J{a,z) J^i Jo {1-xy V [[F{A{T),z';p,P)]^ J{a,z) F{A{T), z'; p, P) 



dx 



2J(a, z) 7-1 

dz 



2J{a,z) 



[F{A{T),z';p,P)]' 
da' 



T ln{l -x)dA{T) 



(1 — xy J{a, z) dx 



T 



^1 - Xi)\E{xi,S,a') 



sgn(£'(xj, S, a') In (1 — Xi) 
J {a, z) 



19 



where we have first integrated the second term of the sum over T by parts (to increase the 
power of F{A{T), z';p, P), the boundary term vanishes when x — 0, 1) and then introduced 
the integration over a' as in the previous section. 

Now, we can use the uniqueness theorem of PTIR [|19| and identify the BS weight function 
on the right-hand side of BSE: 



p^''\a',z') = \ dz daV^^\a',z';a,z)p^^\a,z) 



z'; a, z) 



(A32) 



— 1 — oo 
3 



EE 

T i=± 



d{xi)d{l-Xi)d{D) 



2J{a,zy 



TJ{a, z) 



'I - Xi)\E{xi,S,a 



- - sgn S, a')) ln(l - • (A33) 



Before treating the general case with fuUy dressed propagator we consider in the next sub- 
section the pure ladder BSE. 



a. BSE without propagator dressing (for n = 2) 

In this case we can proceed in a way very similar to the undressed BSE for n = 1, only 
with different re-scaling factor 

p[2](a,^) = J{a,zfp^''\a,z), 

V^'Ka',z\a,z) = 

J [a, zy 

Imposing the auxiliary normalization 

j dz j ap^'^\a,z) = 1, (A34) 



we can write 
P 



(a', z') = X V^J'l {a', z')-X [ dz [ a V}^^ (a', z', a, z) pt'' (a, 2) , (A35) 

^ •' s=± 



i=± [a , z ) 

vp\a',z',a,z) = ^f/; ; j:eix^,)eii-x,)eiD) 

zj [a , z ) 

^ ' ^ sgn S", a')) ln(l — Xj 



Xi=Xi{T=0) 



[1 - Xi)\E{xi,S,a') 



b. Dressed ladder BSE for n = 2 

When the all self-energies are taken into account the function J{a, z) do not factorize 
and numerically convenient redefinition ( |A34| ) can not be used. Nevertheless, we can still 
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separate from the integral over spectral variables and the sum over T the term which 
depends on a and z only through J(q;, z), namely the part for which di — > 5{ai — ml) and 
T = 0. Explicitly, using the notation 

V{a ,z ,a,z,T,ai,a2,a^) = 2^ 

i=± z,J\a^z) 

( TJ{a,z 



—- - sgn {E{xi, S, a')) ln(l - Xi) \ 



and imposing the normalization condition 



1 oo 



-1 -oo ^ ' ^ 

the BSE can be re-written as 

pl'l {a', z') = A Vi^^ {a', z') + X J dz J da^^ V^^^ {a', z', a, z) {a, z) , 

Vi'\a',z') = -l [ dasa{a3)Y.T.^i^iWi^ - ^iWiD)sgn{E{xi,S,a'))ln{l - Xi) , 

^ s=± i=± 

V'^\a' , z' ,a, z) = J da-^dr{a^) J daida2 

[SiS2{l - St,o) + Sia2 + (J1S2 + (71(72] V{a', z', a, z, T, ai, 0:2, a^) , 

where in VQ'^\a',z') we take Xi = Xi{T = 0), ai = ml, a2 = ml and in the last term 
Si — 5{ai — ml) and ctj = (T{ai — m^),i — 1, 2. 
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APPENDIX B: ELASTIC ELECTROMAGNETIC FORM FACTOR 



In this Appendix we derive the expression for the elastic electromagnetic form factor 
The relevant matrix element is diagrammatically depicted in the Fig.^ and the 
starting equation for the e.m. matrix element is given by ( [4 .21 ). The labeling of momenta 
corresponds to the Fig.^: g is a virtual photon incoming momentum and we put = 
—q^ > 0, Pi = P and Pr = P + q are the total momenta of bound state in in- and out-state, 



mo 



respectively. For simplicity, we will consider the equal mass case: mj 

The form factor G{Q^) can easily be extracted from its definition ( |4.1| ). After multiplying 

by Pi + Pf we get 

G,iP + q,P)i2P + qr 



(Bl) 



e(2P + g)2 

The matrix element is evaluated between on-shell states of appropriate composite scalar, i.e 
= (P + g)^ = which imphes 

2Pg + 





The one body current reads 

r{Pf,p^) = {p+2k+qr . 

Using (|B2|) one can simplify 

(2P + qf = + , 
(2P + q) ■ {{2k + q + P) = 2k ■ {2P + q) + 2M^ + 



(B2) 
(B3) 



2 , 



(2P + g) ■ ((2A; + g + P) 
(2P + g)2 

Taking this into account we can write: 

d^k 



1 / 4fc-(2P + g) ' 

2 I 4M2 + Q2 



G{Q') 
HA 

i=l 



(27r)42 



3 

1=1 



rr 



1 + 



Ak ■ (2P + q) 



P P P 

D{k + -; m^)D{-k + -; m^)D{k + ? + "^^) 



rr = T{k + ^,P + q)T{k,P). 



(B4) 

(B5) 
(B6) 



Now, we will express G{Q^) in terms of spectral functions of the bound state vertex 
functions r, re-writing first the product of the propagators with the help of the Feynman 
parametrization. 

For the product of the propagators one gets: 



nA = D{k + q+^)^Jdv 



1 



k'^ + r]k ■ P + ^ — m? + ie 



ds 



dy-, 



2s-l 



[k^ + ^ -m^ - ^Q^ + 2sk- q + yk- P + ie] 



(B7) 
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where the substitution y = s + {1 — s)ri has been used and relations P ■ q = and 
= has been replaced used. 

Next step we combine PTIR (for n = 2) of the product of the bound state vertex functions. 



daida2 / dzidz2P^'^\ai, Zi)p^'^\a2, 



Z2 



dx 



G 



M 



r(4) x(l-x) 

r(2)r(2) 



k^ + — + {zi + x{z2 - zi))k ■ P + x{l + Z2)k ■q-x{l + Zg)-^ 



4 



4 



of 



— {l — x)ai—xa2 + ie. (B8) 

Making use of the Feynman variable t for matching ( [B?] ) with the term in large brackets 
the relation for the form factor (|B4D can be rewritten as 



2r(3) 

d^k 



111 1 

j ds j dyjdxx{l-x) jdtt^{l-t)' 



2s-l 



. I A (2P+q).k 
(2P+g)2 



a 
b 



(27r)4 [c-k^-k-{aP + hq) - lef 
t[zi + x{z2 - zi)] + (1 - t)y , 
tx{l + Z2)t + 2s{l - t) , 

b — — h (1 — t)m^ + t(l — x)ai + txa2 ■ 



(B9) 



where we have omitted the vertex weight functions and the a, z's integrals (exactly the pre- 
factor in front of the large bracket in (p^)). Integration over the momentum k (with the 
shift k + (aP + hq)/2 k then yields 



d'^k 



^ £(7) 

2 r(3) J (27r)4 [c - F - A; ■ (aP + hq) 



I 2k-(2P+q) 



lt\ 



r(5) 

4(4 



i(Pa+qh).(2P+q) 



\2 r 



c+UPa + qby 



(BIO) 



This relation can be further simplified using 

r(4) 



in both the numerator and the denominator 



(4vr) 



111 
dt dx dy 



ds 



x{l-x)t^{l-tf{l-a) 





M2 







-6(6 — a — 1) — r 



(Bll) 



r = — (1— a^)^ — \- {1 — 1)171^ + t(l — x)ai + txa2 , 

a = t[zi + ^(^2 - -2i)] + (1 - , 
6 = tx{l + Z2)t + 2s{l - t) , 

It can be shown that for > the denominator is nonzero. The function I{Q'^) can be 
easily calculated numerically. Including the missing pre-factors the elastic electromagnetic 
form factor is given by: 

roc rl 

G{Q^) = daida2j^dzidz2P^^\ai,zi)p^^\a2,Z2)I{Q^) . (B12) 
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m^/m 


ry = 


r] = 0.2 


r] = 0.5 


T] = 0.8 


T] = 0.999 





1.9998 


1.954 


1.592 


0.9067 


0.03322 


0.5 


2.5663 


2.498 


2.142 


1.421 


0.3873 


Ref.|] 


2.5662 


2.4988 




1.4056 


0.3853 



TABLE 1. Dimensioneless coupling A = g'^/{AnmY as a function of fraction of binding 
rj = \/P^/2m for two cases of exchanged mass m^. The case m^/m = is the 
Wick-Cutkosky modeL The second case m^/m = 0.5 is compared with the result obtained 

by Kusaka et al 0. 





16*16 


32*32 


64*64 


96*96 


CXD 


r] = 0.999; mg = 0.50 


0.3782 


0.3754 


0.3794 


0.3816 


0.3874 


r] = 0.950; mg = 1.00 


1.310 


1.341 


1.355 


1.360 


1.371 


r] = 0.950; mg = 0.50 


0.752 


0.761 


0.777 


0.783 


0.804 


r] = 0.950; mg = 0.10 


0.306 


0.350 


0.375 


0.385 


0.409 


r] = 0.000; mg = 1.00 


3.143 


3.273 


3.342 


3.366 


3.416 


r/ = 0.000; mg = 0.50 


2.207 


2.343 


2.445 


2.483 


2.566 



TABLE 2. The coupling A = g^/lAnmy for bare ladder BSE as a function of the number 

of mesh-points. 



rj 


0.0 


0.4 


0.5 


0.5 


0.5 


0.95 


0.95 




1.0 


0.25 


1.0 


2.0 


4.0 


0.1 


1.0 


A 


3.416 


1.77 


2.928 


4.911 


9.997 


0.409 


1.371 


Ref.|15| 


3.419 


na 


2.940 


na 


na 


0.416 


1.371 



TABLE 3. Dimensioneless coupling A = g'^ / {Airm)'^ for several selections of m^/m 
fraction of binding r] = \/P^/2m, compared to the results of ref. . 



and 
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p+ P/2 




-p+P/2 

FIG. 1: Diagramatical representation of the BSE for the bound state vertex function. 



0.10 



0.08 



0.06 



0.02 - 



0.00 




FIG. 2: The imaginary part of the renormahzed propagators for different values of couphng A, 
calculated from the DSE in bare vertex approximation. Upper curves are for particles $1^2 (which 
have identical self-energies, lower ones for particle $3 .) 




FIG. 3: The dependence of field strength renormalization constants on the coupling A. 
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q 



A" p+q+P/2 p+P/2 



P+q 



\1 



A 



-P+P/2 

FIG. 4: Diagrammatic representation of the e.m. current bound state matrix element. 
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a 

FIG. 6: The re-scaled weight function p{a, z) of the bound-state vertex for 77 = 0.95 calculated in 
bare ladder approximation. The mass of the exchanged boson ms = 0.5m. 




FIG. 7: The same as in previous figure, but for ma = 0.1m and rj = 0.95. 
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FIG. 8: The eigenvalues A calculated for the bare BSE, with dressed kernel V and for dressed 
ladder BSE. Beyond the critical value of coupling Xcrit = 1-5 only the bare solution is available. 
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FIG. 9: The behavior of the elastic electromagnetic form factors for various bound states char- 
acterized by 77. The mass of exchanged particle is fixed to be 1713 = 0.5m, except for the scalar 
deuteron case (D.), which is calculated for comparison using two different grids. 
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■ 2 4 - , 6 8 10 

Q/m 

FIG. 10: Variation of the elastic electromagnetic form factors with the mass of the exchanged 
particle while r) = 0.5 is fixed. The cases rf = 0.4, ms = 0.25m and the scalar deuteron are 
included for comparison. 
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